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Decoupled two-dimensional superconductivity and continuous melting transitions 
in layered superconductors immersed in a parallel magnetic field 
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Computational Materials Science Center, National Institute for Materials Science, Tsukuba 305-0047, Japan 

(Dated: May 14, 2004) 

Possible pliases and the B — T phase diagram of interlayer Josephson vortices induced by a 
magnetic field parallel to the superconducting layers are investigated by Monte Carlo simulations 
based on the anisotropic, frustrated XY model. While for low magnetic fields and small anisotropy 
parameters a single first-order transition is observed similarly to the melting of Abrikosov (or pan- 
cake) vortex lattice, an intermediate phase, characterized by decoupled, two-dimensional (2D) quasi 
long-range crystalline order (QLRCO) and superconductivity, is found at high magnetic fields and 
large anisotropy parameters. Combining the simulation results with a symmetry argument, it is 
revealed that this intermediate phase is of Kosterlitz-Thouless (KT) type, and the melting of 2D 
quasi Josephson vortex lattices and suppression of superconductivity is a KT transition. Evolution 
of the intermediate phase to the low-temperature phase of 3D LRCO is second order and belongs 
to the 3D XY universality class. The three phase boundaries merge at a multicritical point at the 
magnetic field of order B^ic ~ cjf>o/2\/37d^ in the B — T phase diagram. It is revealed that decoupling 
of the 3D Josephson vortex lattice into the 2D phase is triggered by hops of Josephson flux lines 
across superconducting layers activated by thermal fluctuations. The equilibrium phase diagram 
with the KT phase at high magnetic flelds and large anisotropy parameters is consistent with the 
peculiar Lorentz-force-independent dissipation observed in highly anisotropic high- Tc superconduc- 
tor Bi2Sr2CaCu208-i-H by lye et al. (Physica 159C, 433 (1989)). 

PACS numbers: 74.60.Ge, 74.20.De, 74.25Bt, 74.25.Dw 
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I. INTRODUCTION 



The high- Tc superconductors share a layered struc- 
ture in which the superconductivity is widely believed 
to occur mainly in the Cu02 layers intervened by layers 
of charge reservoir (block layers). For many purposes, 
the materials can be considered as intrinsic stacks of su- 
perconducting layers, coupling with each other by the 
Josephson effect. They are therefore strongly anisotropic 
between the normal (c-axis) and in-plane directions. In 
spite of this anisotropy, the superconductors are three- 
dimensional (3D) in nature at temperatures near the su- 
perconductivity transition point at zero magnetic field 
Tc- Especially, the critical phenomena of superconduc- 
tivity transition are governed by the 3D XY universality 
class [l| . Since the coherence lengths diverge in the criti- 
cal region, the short-range (SR) nonuniformity caused by 
the layer structure becomes unimportant. 

There are two cases that physics can be different. 
The first case is realized at low temperatures, at which 
the correlation length ^c becomes comparable with, and 
even smaller than the interlayer separation d, where two- 
dimensional (2D) properties can arise. A crossover be- 
tween 2D and 3D behaviors is expected to occur when 
the value of £,c/d is of the order of unity. The second 
case is realized when the superconductor is immersed in 
a strong magnetic field parallel to the Cu02 layers. The 
coupling between the superconducting layers are weak- 
ened significantly by the magnetic field, and the effective 
correlation length ^c can be very small even at high tem- 
peratures. To elucidate the second case is the objective 
of the present study. As revealed in the present paper. 



there is a genuine thermodynamic phase transition be- 
tween 2D and 3D phases at high magnetic fields. 

It is now well established that the superconductivity 
transition in a type II superconductor in magnetic fields 
[31 is first order, accompanied by the freezing of the flux- 
line liquid into 3D lattice |j,|l,p|. This notion is partic- 
ularly important because the transition was considered 
second order for a long time. We notice however that 
the first-order normal to superconductivity transition is 
mainly observed under magnetic fields perpendicular to 
the layers, where the 2D translation symmetry enjoyed 
by pancake vortices is broken at the transition. In a 
sharp contrast, a parallel magnetic field penetrates the 
layered superconductor in the form of Josephson vortex 
of flux quantum !9| through block layers (2i- The rele- 
vant c-axis translation symmetry is reduced a priori to 
a discrete one, raising the possibility of new phases and 
new melting process. 

As an order competing with the 3D lattice order, de- 
coupled superconducting layers with in-plane quasi-long- 
range (QLR) correlations was proposed by Efetov even 
before the discovery of high-Tc superconductors with pro- 
found layer structures y. Later on, an exponential de- 
pendence of the interlayer shear modulus on the mag- 
netic field is derived by Ivlev et al. ^ , consistently with 
the decoupling scenario since the small shear modulus 
can hardly hold the 3D lattice when thermal fluctua- 
tions set in. This possibility was however questioned by 
Mikheev and Kolomeisky lCr| using the renormalization 
group (RG) approach. It is derived that any weak inter- 
layer coupling should be relevant and only a 3D LR crys- 
talline order (CO) is possible, as far as hops of Joseph- 



son flux lines across Cu02 layers can be neglected, an 
approximation adopted in the previous works. A similar 
conclusion was obtained by Korshunov and Larkin [ll|. 

Transport experiments on high-Tj, superconduc- 
tors have been sheding lights on the possible crys- 
talline order of Josephson flux lines. A peculiar 
Lorentz-force-independent dissipation has been found in 
Bi2Sr2CaCu208+j, under H ~ 5T hy lye et al. that the 
resistivity is independent of the angle between the mag- 
netic field and current when they are both parallel to 
the Cu02 layer [l3|. In the same family of materials, 
non-Ohmic power-law I-V characteristics are observed by 
other groups [ij, [lj| • Chakravarty et al. demonstrated 
that the Lorentz-force-independent dissipation cannot oc- 
cur in a lattice phase |l5| . In order to explain these pe- 
culiar dissipations. Blatter et al. proposed a Kosterlitz- 
Thouless (KT) transition HallS at high magnetic fields, 
relating the creeping of Josephson flux lines across the 
Cu02 layers caused by applied currents with the shear- 
less flux-line state [13 • (See [T^] for a possible KT 
phase at intermediate magnetic fields.) On the other 
hand, an excess resistivity for 7J _L / over H\\I, and 
thus a Lorentz- force- dependent dissipation was observed 
in YBa2Cu307_5 single crystals at low fields |20ll2ill2^. 
Kes et al. addressed this difference in terms of the dif- 
ferent anisotropy parameters |23j |. Taking the Bi-based 
materials almost " magnetically transparent" , the authors 
attributed the experimentally observed dissipations to 
uncontrolled miss alignments of magnetic field. 

Motivated by an experiment by Kwok et al. suggestive 
of continuous melting transition at intermediate mag- 
netic fields 24], Balents and Nelson proposed a smectic 
phase: a regular subset of block layers are selectively oc- 
cupied by Josephson vortices; in the occupied block lay- 
ers Josephson flux lines behave like liquid. They showed 
that melting of the 3D lattice can be continuous in the 
presence of layer pinning [23 . See also |2g . 

Although these studies provide important physical in- 
sights for the mixed states in layered type II supercon- 
ductors, it is clear that discrepancies among different ap- 
proaches have not been resolved and that a unifled pic- 
ture is still not available. The main difficulty lies in the 
features of different energy contributions involved in the 
present system: First, the inter- vortex repulsive force is 
highly anisotropic between the c axis and the ah plane; 
Second, Josephson flux lines feel strong pinning poten- 
tials from the Cu02 layers. At a first glance, approxi- 
mations which take into account differences in strengths 
of these energy contributions can simplify the problem. 
It turns out that the situation is quite complicated. The 
Josephson flux lines arrange themselves such that they 
reside at positions from which deviations in different di- 
rections cost equal energies. Although the profound layer 
structure reduces the c-axis component of the displace- 
ment field at low temperatures, it certainly enhances the 
rigidity of the fiux-line lattice and stabilizes it to high 
temperatures. At the melting temperature, thermal fluc- 
tuations are significant and smear out the strong layer 



pinning to an effectively weak one. Therefore, no ap- 
proximation can be justified easily in the regime of phase 
transition. One thus has to treat the competition among 
the anisotropic repulsion energy, the periodic layer pin- 
ning potential, and thermal fluctuations simultaneously. 

On the other hand, since a Josephson vortex is of large 
size in the in-plane direction, point-like defects play much 
less pinning effect compared with a pancake vortex. It re- 
lieves one from being involved into glassy problems, and 
makes the study on a pure Josephson- vortex system more 
relevant to reality. As a result, the high- field part of the 
B — T phase diagram for a pure vortex system becomes 
accessible experimentally. For this point, please see sev- 
eral interesting proposals pTl |28|, 29] in terms of pan- 
cake vortices, which are however hampered by inevitable 
point-like defects in superconductors which govern the 
behaviors of pancake vortices. 

These situations motivate us to investigate Josephson- 
vortex systems by means of Monte Carlo simulations. 
The Hamiltonian is the so-called anisotropic, frustrated 
XY model on the superconductivity order parameter, 
where the frustrations in phase variables are induced by 
the magnetic field. This model has been used to simulate 
quite successfully the melting phenomenon of Abrikosov 
(or pancake) vortex lattice under magnetic fields paral- 
lel to the c axis jM El El IM HI Isl ■ We adopt the 
same model in the present study on interlayer Josephson 
vortices. The main results are summarized as follows. 
There is a multicritical point in the B — T phase diagram: 
below the multicritical field a single first-order melting 
transition upon temperature sweeping is observed; above 
it there exists an intermediate KT phase characterized 
by in-plane 2D, QLRCO and superconductivity in be- 
tween the normal phase and 3D lattice phase, accompa- 
nied by two continuous melting transitions. The exis- 
tence of a KT phase at magnetic fields above the multi- 
critical point explains the peculiar dissipations, such as 
the Lorentz-force-independent resistance and the power- 
law non-Ohmic I-V characteristics, observed in Bi-based 
materials for which the multicritical magnetic field is 
approximately Smc — 2T presuming 7 = 150. In con- 
trast, Lorentz-force-dependent dissipations are expected 
for YBa2Cu307_5, since to access the KT phase one 
needs a magnetic field above B^ac — SOT taking 7 = 8. 
Some of the results were published in Ref. |3a . W(\ . 

The remaining part of this paper is organized as fol- 
lows. In Sec. II the model Hamiltonian is introduced and 
some details of the simulation techniques are presented. 
After a discussion on the description on flux line lattice, 
we show in Sec. Ill simulation results on the first-order 
melting for a system corresponding to YBa2Cu307_5 sin- 
gle crystals. It then follows numerical evidences on an 
intermediate phase in a highly anisotropic system. De- 
tailed characteristics are provided which allow us to con- 
clude the KT nature of this phase. The superconduc- 
tivity transition is discussed in Sec. IV, and the rela- 
tionship with the crystallin order of Josephson vortices 
is revealed. Section V addresses the possible universality 



of the transition between 2D phase and the 3D lattice. 
The B — T phase diagram of Josephson vortex systems is 
mapped out in Sec. VI. Finally, summary on the main re- 
sults derived from simulations and discussions on recent 
experiments are given in Sec. VII. 



II. MODEL AND SIMULATION DETAILS 

In the presence of magnetic field, the amplitude of su- 
perconductivity order parameter, associated with the lo- 
cal tendency of electron pairing, attains a finite value 
at Hc2 through a crossover. The genuine thermody- 
namic phase transition takes place at a temperature lower 
than that corresponding to Hc2{T). Therefore, the most 
important thermal fluctuations in the thermodynamic 
phase transition come from phase variables of the super- 
conductivity order parameter. Under a magnetic field 
sarallel to the layers, the effective Hamiltonian is thus 



Ti. = —J 2_^ cos{ipi — ipj) 

R.1 — R. j \\x.,y axis 



R-j — R,j II c axis 



R, 



— ^ cos{ipi - ip-j - -;- I A-dR). (1) 



»0 jRi 



The model is defined on the simple cubic grid with the 
unit length equal to the separation between Cu02 layers 
d. The couplings given by J = 0o'^/167r'^A^{, are limited 
to nearest neighboring grid sites and 7 = Xc/Xab- To 
be specific, we put a; _L y _L c and B||y, and choose the 
Landau gauge A = (0, 0, —xB) in Eq.(l). 

The cosine functions of the gauge invariant phase dif- 
ferences in ab planes model the kinetic energy terms in 
a Ginzburg-Landau free energy functional. Higher har- 
monics are neglected since the lowest order terms with 27r 
modulation are sufficient for describing relevant thermal 
fluctuations at phase transitions. The second term in the 
Hamiltonian is the Josephson energy between neighbor- 
ing Cu02 layers. The magnetic field weakens the inter- 
layer coupling by inducing frustrations. Our Hamiltonian 
can be derived from the Lawrence-Doniach free energy 
functional |33|, which is proposed for layered supercon- 
ductors and recognized as a good model for high-Tc su- 
perconductors with profound layer structure. 

Spatial variations and thermal fluctuations in the mag- 
netic induction, and thus magnetic couplings, are ne- 
glected in Hamiltonian (1), which is in principle justi- 
fied only at the strongly type II limit. It is a reasonable 
approximation in the present case. The London penetra- 
tion depth Ac is very large (of order /im), which makes the 
spatial variation of the magnetic induction in the x direc- 
tion and thus the magnetic interactions among Josephson 
flux lines very small. Magnetic couplings may play a role 
when pancake vortices are activated by thermal fluctua- 
tions. It was addressed, however, that the magnetic in- 
teraction between two pancake vortices on a same Cu02 



layer increases with the separation only as IJR at large 
distance because of the interlayer coupling [13 , which is 
therefore much smaller than the elastic energy. There- 
fore, magnetic couplings can be ignored safely as far as 
the thermodynamic phase transition is concerned. 

Thermal fluctuations in the amplitude of superconduc- 
tivity order parameter can be taken into account in sev- 
eral ways. One can include a temperature dependence of 
the amplitude in the mean-field way, such that the pene- 
tration depth and the coupling constants vary with tem- 
perature. This treatment permits a quantitative compar- 
ison between simulated quantities and experimental ob- 
servations [32, I33 • Alternatively, one can take thermal 
averages of superconductivity order parameter on scales 
larger than the unit length of grid in ab planes. A fully 
disordered configuration of phase variables on the scale of 
d then results in zero superconductivity order parameter 
at any larger length scales. 

The Hamiltonian (1) provides a reasonable description 
for the layer pinning on Josephson vortices in the whole 
temperature regime, even without taking into account 
the temperature dependence of Ginzburg-Landau corre- 
lation length. At low temperatures, the layers modeled 
by the discreteness of the underlying grid in the c direc- 
tion set a series of barriers to the motion of Josephon 
flux lines. If temperature is high enough, on the other 
hand, pancake vortices are activated thermally on Cu02 
layers, such that Josephson flux lines can hop to neigh- 
boring block layers easily. Thermal fluctuations reduce 
effectively the layer pinning potential. 

The above Hamiltonian has a trivial limit 7 = cx), 
where KT transition takes place at T^^"^ ~ O.SQJ/ks 
in each independent ab plane. For large but finite 
anisotropy parameters, there is a broad crossover regime 
which is dominated by this limit, which can be seen in 
the specific heat |37ll40l |. 

In the present study, the magnetic field is fixed at 
/ = Bd^ /(j)o = 1/32 while the anisotropy parameter is 
tuned for the convenience of simulation. The system 
size is Lx X Ly X Lc = 384d x 200d x 20d under peri- 
odic boundary conditions. There are totally Ny = 48000 
Josephson vortices in the ground state. The system size 
is set anisotropically since we are interested in systems 
of large anisotropy parameters: there are / x L^/d = 12 
Josephson flux lines the x direction, which is comparable 
with the number of Josephson flux lines in the c axis: 
Lc/2d = 10. (See Fig. 2 oi '37^.) As the field direction 
coincides with ab planes of strong couplings, the system 
size in this direction has to be taken large enough in or- 
der to treat thermal fluctuations sufficiently. Although a 
full analysis on finite-size effects is still difficult with the 
available computing resource, it will be revealed later by 
comparing simulation results for different anisotropy pa- 
rameters that the phases and phase diagram derived in 
the present paper are free of serious finite-size effects. 

In a typical simulation process, we start with a random 
configuration of phase variables at a sufficient high tem- 
perature. The system is then cooled down gradually with 



the temperature skips and the number of Monte Carlo 
sweeps hsted in Ref . |37l| . After arriving at the lattice 
phase at a sufficiently low temperature, T = OAJ/ks, 
we heat the system back slowly. Thermal averages are 
taken over ~ fO^ Monte Carlo sweeps at temperatures 
around the transition points. All the results shown in 
the present paper are calculated in the heating process. 



III. CRYSTALLIZATION TRANSITION 

A. Description of Josephson vortex lattice 

The crystalline order of the system is described by the 
correlation function 



5(R) = (p(R)p(0))-(p(R))(p(0)), 



(2) 



where R = (x, y, z) and /3(R) is the y component of the 
vorticity at position R which is explicitely defined by 
/9(R) = J2a ^(I^-L ~ ^±,a{y)) with Rj_ — {x, z) , and its 
Fourier transformation, i.e. the structure factor. 



The cross terms between different reciprocal vectors dis- 
appear because of the infinite expectation values of dis- 
placement fields. The structure factor is then given by 



S{k) ^ J2 <i'R[e~'^''-'^^y^ + e-^(''+K.)-R]5(K,,R). 
j •' 

(8) 

For the present flux line system, it is convenient to 
consider partial structure factors in the section per- 
pendicular to the magnetic field. They are related to 
the above 3D one by a partial integral, S{kx,y,kc) ~ 
j dkye^^y"^ S{}<i). In the present paper, we will show sim- 
ulation results on S{kx,kc) = S{kx,y = 0, fcc) which re- 
veals the vortex correlations in the same crosssection per- 
pendicular to the magnetic field. 

In a 3D crystal phase, (7(Kj,R)'s approach to con- 
stants at large distances, and (5-function Bragg peaks 
should be observed at the reciprocal-lattice vectors. For 
the present system of weak interlayer couplings, there is a 
competing order to the 3D crystalline order, in which the 
Josephson vortices are strongly correlated in each block 
layer such that a QLRO is present, while only SR corre- 
lations are realized in the c direction: 



5(k) 



d^Re 



-ik-R 



5(R). 



(3) 



In lattice phase subject to thermal fluctuations, the 
vortex density is expressed in terms of the reciprocal vec- 
tors as 



p(R) = *oe 



-iKo R 



-*i(R)e 



jKi R 



-*2(R)e"'^^-^+c.: 



(4) 
where K^ with 7 = 0, 1, 2 arc primitive reciprocal lattice 
vectors and Kp-I-Ki -I-K2 = 0. Higher harmonics are not 
included for simplicity, which are not important for the 
later discussions on phase transitions. ^j(R)'s are order 
parameters of the crystalline order, and are expressed by 
the displacement field u = (u^jjO, Uc) pdl |: 



*j(R) = e^^r^i^). 



(5) 



The density-density correlation is then given by corre- 
lation functions of order parameters 



sm =Y.M 



iKjR, 



^,{K)^m)+c.c:_ 



with the correlator of displacement fields defined by 



(6) 



7(K,R) 



-iK-[u(R)-u(0)] 



)• 



(7) 



5(K,,R)^e 



-\A/ic 



/r" 



(9) 



with R — (r, z), and r an in-plane positional vec- 
tor rescaled from the original one according to the in- 
plane elastic constants [l^ |2a|- As will be revealed 
later, this case occurs when the mangeic field and the 
anisotropy parameter are large enough, such that every 
block layer is occupied by Josephson vortices and the 
two unit vectors in real space are ai — {d/ f, 0, 0) and 
3-2 = (f^/S/, 0, d). The primitive reciprocal lattice vec- 
tors are K = (±2/7r/d, 0,±7r/d) and (0,0, ±27r/d). The 
last two ones are equivalent to the origin. In this case, 
dominant thermal fiuctuations are associated with the 
other four reciprocal vectors. Especially, the fee and k^ 
profiles of the partial structure factor with y = around 
them are evaluated as 



5(±2/7r/d, kc) - / da;dze-^(*=<^±^/'"^5(K, R) 



C^ 



i--' + {k,±7r/d)^ 



^\kx±2fTr/d\-^+^, 



(10) 



(11) 



with R — (a;,0,z), up to multiplicative coefficients and 
linear corrections. Namely, in such a phase of 2D QLRO 
the structure factor should show power-law singularities 
in the k^ direction at reciprocal-lattice vectors, while 
smears out in the k^ direction in a Lorentzian form. By 
fitting the profiles of the Bragg spots, the correlation 



length ^c and the exponent rj of QLRO in ab planes can 
be evaluated. Detailed discussions on profiles of Bragg 
peaks in the smectic phase of liquid crystals can be found 
inRef.El. 



B. First-order melting transition: 7 = 8 

Let us start with a system of the anisotropy param- 
eter 7 = 8, which models the high-T^ superconductor 
YBa2Cu307_5. Upon temperature sweeping, a first- 
order transition is observed, indicated by the (5-function 
peak in the specific heat and a kink anomaly in the 
Josephson energy at r,„ = 0.96J/fcs, as displayed in 
Figs.l and4of Ref. 37]. 

In order to characterize the phase transition better, 
we investigate the crystalline order of Josephson vortices 
in terms of the structure factors. At low temperatures, 
there is a fine 3D LRCO in the system indicated by the 
(5-function Bragg peaks at the reciprocal-lattice vectors, 
not only the minimal ones but also at higher-order satel- 
lites, as displayed in the top-left panel of Fig^ The 
Bragg peaks are observed up to T^ — 0.96 J/ /cs, without 
losing their sharpness as shown in the top-right panel of 
FigQand in FigO The 3D LRCO disappears when tem- 
perature crosses Tm- The first-order transition at T„j is 
therefore identified as the melting of the Josephson vortex 
lattice, similar to the melting of Abrikosov (or pancake) 
vortex lattice in magnetic fields parallel to the c axis. 

At T ~ 0.98J/fcsi SR correlations are observed as 
shown in the bottom-left panel of Fig^ which charac- 
terize a vortex liquid with layer modulation. It is wor- 
thy to notice that, in spite of the high anisotropy in the 
vortex lattice, thermal fluctuations of large wave-lengths 
are quite isotropic for 7 = 8 as reflected in the spots at 
{kx,kc) = (±2/7r/(i, ±7r/(i) in the bottom- left panel of 
FigH At higher temperatures, T > 1.2J//cb, correla- 
tions are suppressed by thermal fluctuations to the scale 
of grid size of simulation as depicted in the bottom-right 
panel of Fig^ Containing huge number of closed vor- 
tex loops the system is completely disordered. In such 
a case, a coarse-graining in the ab planes results in zero 
amplitude of the superconductivity order parameter, cor- 
responding to the normal states. In this temperature 
regime vortices are therefore not suitable for describing 
the system anymore. The transformation from the vortex 
liquid to the normal phase is a crossover, characterized by 
the broad peak in the specific heat in Fig. 3 of Ref. |37| . 

We have also simulated systems with 1 < 7 < 8 at 
/ = 1/32. The melting point increases with decreas- 
ing anisotropy parameter. Josephson vortices distribute 
in a subset of block layers for 7 < 5 [33. The phase 
transitions are first order at which 3D LRCO sets up for 
all the cases we have searched. These simulation results 
look inconsistent with scenarios of continuous melting at 
intermediate and low magnetic fields. 




FIG. 1: Structure factors S{kx,kc) for 7 = 8 at several 
typical temperatures. Top-left: T — O.lJ/fcs; Top-right: 
T = 0. 96 J/fcs; Bottom-left: T = 0. 98 J/fcs; Bottom-right: 
T = 1.5J/fc_B. The panels are for wave numbers within k^ G 
[-257r/192d, 257r/192d] (horizontal) and kc G [-27r/d, 2-K/d] 
(vertical). The spots at {kx,kc) — {0,±2n/d) are quivalent 
to (0, 0) in the present system. 
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FIG. 2: kx (left) and kc (right) profiles of Bragg peaks in the 
top-right panel of FigQfor 7 = 8 at T = 0.96J/fcs. 



C. Multicritical point 

The first-order melting of Josephson vortex lattice is 
suppressed to continuous transitions when the anisotropy 
parameter increases to 7 = 10, when the filling factor is 
fixed at / = 1/32, as reported in Ref. |33- A critical 
value of the product between the anisotropy parameter 
and filling factor may be given by jl^ ^T^ : 



fl 



2V3' 



(12) 



It is easy to see that below the critical value, the 
Josephson vortex system is physically equivalent to a 3D 



anisotropic Abrikosov vortex lattice. Rescaling the sys- 
tem by the anisotropy parameter should result in an ef- 
fective isotropic system. All physical properties, includ- 
ing possible phases and the nature of phase transition, 
should be essentially the same as the isotropic system 
except for an anisotropy scaling |43|. This physical dis- 
cussion is consistent with our observations presented in 
the preceding subsection. 

On the other hand, the system becomes un-rescalable 
above the critical value, since the minimal inter-vortex 
separation in the c axis 2d has been reached at the crit- 
ical value. Any increase in either / or 7 results in over- 
whelming x-direction inter-vortex repulsions, which en- 
hance the in-plane crystalline order. Although the Cu02 
layers fix the inter- vortex distance for Josephson vortices 
in the c axis, it has nothing to do with the inter- vortex 
correlations. This is the most important peculiarity of 
the interlayer Josephson vortices compared with liquid 
crystals in which smectic orders are realized. 




D. Intermediate phase: 7 = 20 

In order to reveal the crystalline order of Josephson 
vortices above the critical point, we choose 7 = 20 which 
is quite above the multicritical value 7 = 16/\/3 ~ 9.24 
at / = 1/32. An appropriately large anisotropy param- 
eter realizes the intermediate phase in a wide tempera- 
ture regime while leaves a reasonable interlayer Joseph- 
son coupling, and thus is convenient for computer sim- 
ulations. For T < 0.65J/kB, a 3D LRCO same as that 
for 7 = 8 is observed in this highly anisotropic case as 
shown in the top-left panel of FigEland in Fig01 

In a sharp contrast to the case of 7 = 8, however, the 
Bragg spots for 7 = 20 at {kx,kc) — (±2/7r/d, ±7r/d) 
smear significantly in the kc direction as temperature in- 
creases to T = 0.7J/kB, while leaving the sharpness in 
the k,j: direction almost unchanged as seen in the top- 
right panel of Fig|31 Thermal fiuctuations are therefore 
very anisotropic as well as the structure of the vortex 
lattice. The diffusive and stripe-like Bragg spots survive 
to T ~ 0.95 J /ks (see the bottom- left panel of FigEJ. 

As shown in the right panel of Fig|Sl the kc profile 
of the Bragg spots at T = 0.7 J/kg is fitted well by 
the Lorentzian function in Eq. ljlOII with the correlation 
length ^c — l-5(i. Therefore, the Josephson vortices are 
coupled to each other only SR in the c axis. Although 
it is not easy to distinguish numerically the kx profile in 
the left panel of FigElfrom a (5-function Bragg peak, it is 
easy to see that a (5-function one is impossible as far as 
the system is decoupled in the c axis. 

The kx profiles of structure factors S{kx, kc) for higher 
temperatures are depicted in Figl^l They arc clearly dif- 
ferent from the (5-function peak for T < Q.GbJ/ks for 
the present system, and those for 7 = 8 in the whole 
temperature regime T < Tm- Supposing a cusp singular- 
ity in the profiles in Fig|Sl we can evaluate the exponent 
rj in the correlation function Eq.(|5Jl by fitting the data 



FIG. 3: Structure factors S{kx,kc) for 7 = 20 at several 
typical temperatures. Top-left: T = 0. 65 J/fcs; Top-right: 
T = 0.7.7/ ks; Bottom-left: T = 0.95 J/fcs; Bottom-right: 
T — 0.97 J/fcs. The panels are for wave numbers within kx G 
[-257r/192d, 257r/192d] (horizontal) and fcc G h^n/d, 2-K/d] 
(vertical). The spots at {kx,kc) — {0,±2n/d) are quivalent 
to (0, 0) in the present system. 
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FIG. 4: kx (left) and fcc (right) profiles of Bragg peaks in the 
top-left panel of FigElfor 7 = 20 at T = 0.65 J/fcs- 
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FIG. 5: kx (left) and fcc (right) profiles of Bragg peaks in the 
top-right panel of Fig^lfor 7 = 20 at T = 0.7J/fcs. 
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FIG. 6; kx profiles of structrure factors S{kx, kc) around tlie 
Bragg spots for 7 = 20 in tfie intermediate pliase. Tfie solid 
curves are results of the least-squares fittings to the power- law 
function as described in text. 



points including the peak values. This analysis results 
in 7? ~ 1.97 ± 0.07, 1.42 ± 0.01, 1.17 ± 0.007, 1.11 ± 0.003 
at T = 0.96, 0.95, 0.92, 0.8 J/fcs, with the error bars from 
the least-squares fittings. The fitting curves are shown 
in FigElby the solid curves. The good quality of the fit- 
tings with T] > \ for all the temperatures, especially for 
^ < kx < 7r/16d, can be taken as a justification of the 
assumption on the cusp singularity. 

The temperature dependence of the maximal value of 
Bragg peaks is shown in Fig[7| The intensity of Bragg 
spots is small and increases very slowly with decreas- 
ing temperature for 0.7J/fc_B < T < 0.95J//cb. Across 
Tx where the spots become (5-function like, the intensity 
increases noticably. As a comparison, that for 7 = 8 in- 
creases very sharply with decreasing temperature as soon 
as they become observable at the melting point Tm as 
shown in Fig[7| 

One can also evaluate the partial structure factor 
S{kx,ky) = S{kx,ky, z = 0), which describes the Joseph- 
son vortex correlations in the same block layers. The 
profiles around the Bragg spots {kx,ky) = (±2/7r/d, 0) 
should behave as 



S{kx,0)^\kx±2f7r/d\-^+\ S{2f7r/d,ky)^\ky\- 



2+v 

(13) 
The simulated results are displayed in Fig|Sl for T = 
O.SJ/ks- A power-law singularity is observed, which is 
consistent with that in the x — c crosssections. The col- 
lapse of two profiles left to the peak in Fig|Sl indicates 
that the in-plane correlation functions are governed by 
the same exponent rj. The deviation on the right side 
of the peak is obviously caused by the second peak in 
the kx profile at kx = 4/7r/(i, which is absent in the ky 



FIG. 7: Temperature dependence of the intensity of Bragg 
peaks for 7 = 8 and 20. 
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FIG. 8: Profiles of the structure factor S{kx, ky) around the 
Bragg spot ikx,ky) = (7r/16d,0) for 7 = 20 at T = 0.8J/fcs. 



direction. Since 77 < 2 observed in the structure factor 
S{kx, kc), the unnormalized profiles should diverge in the 
thermodynamic limit as in the above equations. There- 
fore, one cannot evaluate the exponent 77 directly from 
the profiles. 

The above observations therefore indicate the existence 
of an intermediate phase for 7 = 20, characterized by 
SR interlayer correlation and in-plane QLRO, which is 
clearly absent for 7 = 8 below the multicritical point. 

Is an intermediate phase plausible in the present sys- 
tem theoretically? According to the perturbative RG 
expansion to a 2D system by Balents and Nelson |25l |. 
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FIG. 9: Temperature dependence of ratios of Josephson flux 
lines which contain segments hopping into neighboring block 
layers, of those which either hop or collide with neighbors in 
the same block layers, and of populations of thermally excited, 
closed loops of Josephon vortices (normalized by 40 x 240) and 
of those containing also pancake vortices (normalized by 240). 



the lattice is stable at low temperatures with ry < 1/2 
where dislocations caused by hops of segments of Joseph- 
son flux lines across the superconducting layers are irrel- 
evant while the interlayer coupling is relevant; at high 
temperatures with rj > 2 where the dislocations are rel- 
evant while the interlayer coupling is irrelevant, the sys- 
tem behaves as liquid. In the intermediate temperature 
regime, both of them are relevant and the fate of the sys- 
tem is not very well controlled. The system may take one 
of the two states, or it can take the smectic order as the 
authors suggested. 

In simulations, we analyze the trajectories of the 
Josephson flux lines when they travel through the sam- 
ple. As depicted in FigEl Josephson flux-lines are com- 
pletely confined by the Cu02 layers below Tx . Without 
dislocations caused by the hops, the system is in lat- 
tice phase as shown in the top-left panel of Fig|31 consis- 
tent with the theory. Dislocations become popular in the 
system when temperature is elavated cross Tx as shown 
in FigO It is interesting to observe that the maximal 
value of the exponent rj in our simulations, 77 ~ 1.97 
at T ~ 0.96J//c_B, is close to the theoretical prediction 
?7c = 2 for the phase boundary between liquid and the in- 
termediate phase. However, we should notice that while 
the theory predicted relevance of interlayer coupling in 
the intermediate regime, our simulations indicate the SR 
c-axis correlation. The origin of this discrepancy is not 
clear at this moment. As pointed out in Ref.[i3, one pos- 
sibility is that although the long-wave-length displace- 
ments of Josephson flux lines in the neighboring layers 
are coupled, accumulation of disorders in a stack of lay- 
ers leads to decoupling, namely SR c-axis correlation. 



In Fig|2| we also display the temperature dependence 
of population of closed Josephson vortex loops activated 
thermally. These excitations have not been treated in the 
elastic theory, and are responsible for the temperature 
dependence of the exponents rj. The elastic constants in 
the present system are therefore not simply dominated by 
entropy as discussed in the elastic theory. Our system is 
similar to the thermotropic smectics in this aspect [42] . A 
counter example is found in the lyotropic lamellar phase 
of a multilayer fluid membrane system [4J| . See ^^ for 
a detailed discussion. 



E. Novel Kosterlitz-Thouless phase 

Since the density-density correlations are SR in the 
c direction in the intermediate phase, the effective free 
energy is governed by the long-wave-length fluctuations 
given in the following density expression 

p{r) = ^/j{r)e-"^°-'' + V'*(r)e*•'''■^ (14) 

with ko = (kxo, kyo) = (2/7r/ci, 0) and 



V;(r) = IV'le*^-^^"-/^ 



(15) 



When a small amplitude |V'| is set up [23, presumably 
at T ~ 0.96 J/fcs as in FigsEl and [7| the effective free 
energy is reduced to 



(16) 



F = / d'r{Vu,y 



Since the displacement field Ux is continuous and of the 
modulus d/ f, the above free energy is effectively the same 
as the Hamiltonian of the 2D XY model. The intermedi- 
ate phase therefore falls onto the KT fixed line. 



IV. SUPERCONDUCTIVITY TRANSITION 

So far we have concentrated on the crystalline order 
of Josephson vortices. It is also important to investigate 
the superconductivity in the system. Although the am- 
plitude of the local order parameter of superconductivity 
is fixed in the Hamiltonian (1), the fate of macroscopic 
superconductivity is governed by the order of Josephson 
vortices. We measure the LR superconductivity by the 
helicity modulus ^h, ^] . 

As shown in Fie llOL in-plane helicity moduli for 7 = 8 
set up sharply at the melting point T^ ~ 0.96J/A:b, 
where 3D LRCO is realized. The finite helicity modu- 
lus Ty along the magnetic field indicates the breaking of 
U(l) gauge symmetry, and thus the appearance of LRO 
of superconductivity. The finite Ta, is, on the other hand, 
the signature of the intrinsic pinning on Josephson vor- 
tices from the layer structure, since the helicity modulus 
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FIG. 10: Temperature dependence of in-plane helicity mod- 
uli for 7 = 8 and 7 = 20. The solid line is for theoretical 
prediction of the universal jump at KT transitions. For a 
pure 2D unfrustrated XY model corresponding to 7 = 00 in 
Hamiltonian (1) T^?^" ~ 0.89 J/ kg. 



Tx measures the energy cost for sliding Josephson vor- 
tices in the c direction. Tc is vanishing down to zero tem- 
perature, reflecting the absence of pinning force in the x 
direction. Therefore, the system is superconducting in ab 
plane for T < Tm- The sharp drop of the in-plane helic- 
ity moduli at the melting point is in accordance with the 
first-order nature of the melting transition. 

A typical configuration of phase variables on an ab 
plane at low temperatures is displayed in Fig llll The 
phases are uniform along the magnetic field. It is also reg- 
ular in the x direction, except for a single-wave-number 
modulation governed by the density of Josephson flux 
lines lillil. 

The static phase wave can be displayed in a more trans- 
parent way by showing the supercurrent in the x direction 
/a;(r) ~ sin((/3(r -I- x) — f{Y)). The supercurrents associ- 
ated with the phase distribution in Fig llll are displayed in 
the left panel of Fig^| Regions of positive and negative 
supercurrents in the x direction are decorated by white 
and black tones respectively. There are four black (and 
white) stripes in the region of l^ = 128d, corresponding 
to the density of Josephson fiux lines / = 1/32. The pat- 
tern of supercurrents on a neighboring ab plane is shown 
in the right panel of Fig^l The stripes on the two planes 
are opposite to each other in black and white tones, and 
sustain four Josephson flux lines between them shown in 
Fig.2of Ref.0. 

Variations of the phase variables along the c direction 
is displayed in Fig^| The configuration is similar to 
the sequence (0, 0, tt, tt, 0, 0, • • •) derived in Ref.|2g. The 
simulated c-axis modulation is more complex, which to- 







FIG. 11: A typical configuration of phase variables on an 
ah plane at T = O.lJ/kB for 7 = 8. Shown is a region of 



dimensions Ix x /„ 



: 64d X 64d. 




FIG. 12: Supercurrents between two sites neighboring to 
each other in the x direction. Shown is the region of the same 
origin of Fig llll but of larger dimensions l^xly = 12Sd x 128d. 



gether with that in the x direction minimizes the energy 
for weak but finite inter layer couplings. 

The amplitude of the phase modulation in x direction 
is small around and above the multicritical point in accor- 
dance with the analytic result by Bulaevsky and Clem: 
A = l/27r(7/)2 < 1 [i^. Most of the frustrations in- 
duced by the magnetic field are confined in the c direction 
as clearly seen in Figllll The large phase modulations 
in the c direction make the phase slip hard, resulting in 
the finite helicity modulus T^,. The small, and thus soft 
phase modulations in the x direction suppress the helicity 
modulus Tc to zero. 

For 7 = 20, finite in-plane helicity moduli signal the 
onset of QLRO of superconductivity, as a twin of the QL- 
RCO of Josephson vortices. The sharp drops of helicity 
moduli for 7 = 20 correspond to the universal jump in 
KT transitions 1^. 

The jump of helicity modulus in Fig llOl is rounded 
from the universal one by the finite-size effect as usual. 
An analysis on finite-size effects is certainly helpful here. 
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FIG. 13: A typical configuration of phase variables on an ac 
section at low temperatures for 7 = 8. Shown is a region of 
dimensions Ix x Iz ~ 64d x 20d. 



other. It is reasonable since, first, both first-order and 
KT transitions reveal sharp drops in the helicity mod- 
ulus, and secondly, the anisotropy parameters are both 
very large such that the transition points are close to 
the lower bound T^?^'' = 0.89 J/ /ts. The hehcity mod- 
uli in the x and y directions are close to each other for 
7 = 8 merely because that it is close to the multicritical 
point. For smaller anisotropy parameters, thus far from 
the multicritical point, we find clearly that the isotropy 
in the in-plane helicity modulus is broken. This explains 
the anisotropic resistivity in samples of YBa2Cu307_5. 



However, as one can see in a careful study by Lee and 
Teitel on a 2D Coulomb gas system in Ref . |48l| . this ap- 
proach is usually very hard since the finite-size effects 
are logarithmically weak [42j. A similar analysis for the 
present system should be even harder, since first, there is 
a SR correlation in the c direction, secondly, the system 
size should be very huge in order to contain sufficient 
number of Josephson vortices, and thirdly, the present 
system is anisotropic in all three directions. 

For these reasons, we will not seek to verify the uni- 
versal jump in the present system. Instead, upon reveal- 
ing the existence of the novel KT phase from simulation 
results and the symmetry argument in terms of the crys- 
talline order of Josephson vortices in the preceding sec- 
tion, we adopt the expected universal jump [43| 



T/ksTKT = 2/7r. (17) 

to estimate the KT transition temperature l50i |. From 
data in FigUHl we find Tkt ^ O.QlJ/fcs. 

The KT transition temperature for 7 = 20 is above 
that of the isolated 2D system T^^" ~ 0.89J/fcs. This 
is physical since couplings between nearest neighboring 
layers enhance ordering. The bare KT transition point 
sets a lower bound for the normal to superconductivity 
transition in magnetic fields parallel to the ab plane, as 
can be read from Hamiltonian (1). 

An important feature in the helicity moduli for 7 = 20 
is revealed in Fig^| Namely, the helicity moduli in x 
and y directions are collapsing at all temperatures. Pre- 
suming that the relation between the exponent a{T) de- 
fined in the I-V characteristics V ~ /"^t^) and the helicity 
moduH IMIEi 



1 + TpTr/fesr, for fi = x,y 



(18) 



in the KT phase derived in pure 2D systems is also ap- 
plicable for the present KT phase, the collapse of he- 
licity moduli in the two in-plane directions is consis- 
tent with the orientation-independent, i.e. Lorentz-force- 
independcnt, dissipations when the current and magnetic 
field are both parallel to the ab plane [IJ . Further inves- 
tigation is expected in order to address this point com- 
pletely. 

One might notice that in Fig^| the helicity moduli 
for the two anisotropy parameters are very close to each 



V. PHASE TRANSITION AT T^ 

What is the nature of the phase transition between the 
2D QLRO and 3D LRO at Tx ? In order to answer this 
question, we follow Balents and Nelson to compose an ef- 
fective Landau free-energy functional for the 3D Joseph- 
son vortex lattice, which was formulated for a possible 
smectic to crystal transition J23. The vortex density 
in the 3D crystal phase is expressed by three primitive 
reciprocal-lattice vectors as in Eq. Q . A Landau free en- 
ergy for the dominant fluctuations is given by, up to a 
spatial anisotropy in the coefficients which is unimpor- 
tant here. 



+3(V*i V*2 + V*i V*2) 



^(1*1 



|*2P)+W(*1*2 + 1'1*5) 



(19) 

where the order parameter \fo = 1 associated with 
{kx,kc) — (0,±27r/(i) has been included into the coef- 
ficients g and w. Our simulations have revealed that 
Josephson flux lines are completely confined by super- 
conducting layers, namely Uc = 0, for T <T^. Thus, the 
two order parameters in Eq.(|SJl become complex conju- 
gate to each other: 



The free energy is then reduced to 



F 



d?R{^u^f 



(20) 



(21) 



The phase transition is therefore second order with the 
critical phenomena governed by the universality class of 
the 3D XY model. 

A phase transition in this universality class possesses 
a negative critical exponent a for the specific heat. This 
naturally explains the smooth temperature dependence 
of the simulated specific heat around Tx in Fig. 3 of 
Ref.il. 
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FIG. 14: B — T phase diagram for inter layer Josephson vor- 
tices with a multicritical point at Bmc = 0o/2-\/37d'^- The 
phase boundaries Tm{B),TKT{B) and Tx{B) are associated 
with first-order, KT and 3D XY phase transitions as discussed 
in text. Possible crossovers are included by dashed lines. 



When the critical point Tx is approached from below, 
interlayer shear modulus of the Josephson vortex lattice 
is softened continuously to zero as Cqq ~ d/^c ~ {Tx — 
TY with v ~ 0.67, and ramains vanishing in the whole 
intermediate phase. 



VI. PHASE DIAGRAM 

Based on the analyses presented so long, we map out in 
Fig^^the B ~T phase diagram for interlayer Josephson 
vortices, noting that the same physics should occur when 
the magnetic field is tuned at a given anisotropy param- 
eter. The first-order melting line at low magnetic fields 
branches into two phase boundaries at the multicritical 
point characterized by the magnetic field 



^tll.C. 



2^37^2 



(22) 



The phase at high magnetic fields and intermediate tem- 
peratures is a novel KT phase. 

There are two dashed, crossover curves in the phase di- 
agram Fig^] The high-temperature one denoted by Tc2 
follows the broad cusp in the specific heat where huge 
numbers of vortex loops are excited. It is a crossover 
where the amplitude of superconductivity order parame- 
ter attains a finite value, and thus corresponds to Hc2 in 
the mean-field theory. The low-temperature and almost- 
field-independent line Th denotes a crossover tempera- 
ture below which thermally activated pancake vortices 
are very few and bound to each other too tightly such 
that Josephson flux lines cannot hops across the Cu02 
layers. Transport properties can be different on the two 
sides of the crossover line Th since randomly distributed 



point-like defects might pin Josephson flux lines via the 
pancake vortices present above Th |5a. Im ) . As the dis- 
placement in the c direction is massive below Th, one 
expects the transverse Meissner effect in the low temper- 
ature regime. 



VII. SUMMARY AND DISCUSSIONS 

The main finding of the present work is that the normal 
to superconductivity phase transition in magnetic fields 
parallel to superconducting layers is first order at low 
fields, while two-step and continuous with a KT type 
intermediate phase at high fields. There is a multicritical 
point at magnetic field of order B^c — 't'o/'^vild? in the 
B — T phase diagram. 

That the above phases are genuine for the interlayer 
Josephson vortices rather than artifacts from finite-size 
and short-time effects can be concluded from the follow- 
ing considerations: First, the successful observation of 
the single, first-order melting transition for 7 = 8 indi- 
cates that the system size is sufficient for the onset of 
3D LRCO below the multicritical anisotropy parameter; 
Second, the SR order in the c direction for 7 = 20 is not 
a finite-size effect, since periodic boundary conditions are 
adopted in the present study, which tend to enhance or- 
dering. It is not a short simulation-time effect either, 
since a heating process from a 3D LR crystalline order 
is adopted; Third, it is clear that the system described 
by Hamiltonian (1) with a finite 7 should order (maybe 
partially) at a temperature higher than that of the limit 
case 7 = 00, which composes of independent 2D super- 
conductors with T^° ~ 0.89J/kB- Thus, the in-plane 
vortex correlation functions at, for example, T = 0.7 and 
0.8J/kB < T]^T° for 7 = 20 cannot be SR; with the 
short rangeness of c-axis correlations, they cannot be LR 
either. They then should be QLR, consistent with that 
demonstrated numerically. 

The peculiar transport properties observed in 
Bi2Sr2CaCu208+y when magnetic fields and currents 
are both parallel to the Cu02 layers are able to be 
explained in terms of the B — T phase diagram, since 
dissipation at small currents should be governed by 
the equilibrium properties. In the KT phase, isotropic 
helicity moduli with the universal jump at the transition 
point are responsible to the isotropic, and therefore 
Lorentz-force-independent, dissipation and the power-law 
non-Ohmic I-V characteristics. In the 3D lattice phase, 
which is realized at low magnetic fields for Bi-based 
materials or up to quite strong fields in YBCO materi- 
als, larger dissipations should be observed at transverse 
currents than parallel ones. 

The B — T phase diagram in Fig^^ is consistent 
with several recent experiments. By detecting the os- 
cillation in the c-axis flux flow of Josephson vortices in 
Bi2Sr2CaCu2084_j^, Ooi and Hirata succeeded in mea- 
suring the phase boundary on which the 3D triangular 
lattice softens 55]. The resultant curve has the same 
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shape formed by Ty{B) at high fields and Tm{B) at 
low fields, with a kink corresponding to the multicriti- 
cal point. In the same material, Mirkovic et al. observed 
that the sharp drop of the resistivity associated with the 
first-order vortex lattice melting at low magnetic fields is 
suppressed into a smooth one when the magnetic field is 
elevated to about IT, suggesting a continuous phase tran- 
sition 56] . The steep normal to superconductivity phase 
boundary at high magnetic fields observed by Lundqvist 
et al. 1^ is in accordance with the lower bound T^^ on 
Tkt{B), although the KT features are still not available 
experimentally. 

Recently Kakeya et al. observed two plasma modes 
when applying a magnetic field parallel to the ab plane, 
of frequency higher (lower) than that of zero magnetic 
field and increasing (decreasing) with the magnetic field 
[58|. While the high branch is caused by the resonance 
between the electromagnetic field and coherent motion of 
interlayer Josephson vortices, the low branch is asigned 
to the shear vibration of Josephson vortex lattice |59l| . 
Since the KT phase proposed in the present work is char- 
acterized by vanishing interlayer shear modulus, the low 
plasma mode should disappear and the high mode is still 
observable as the magnetic field or temperature approach 
the phase boundary Tx from the 3D lattice phase. There- 
fore, the Josephson plasma phenomenon provides a pow- 
erful technique to test our phase diagram. 

Since the multicritical field is approximately SOT for 
YBa2Cu307_5 presuming 7 = 8 and d = 12A, the nor- 
mal to superconductivity phase transition should be first 
order at magnetic fields available in laboratory to date. 
This is consistent with a recent measurement on the spe- 
cific heat by Schilling et al. up to lOT [s3| ■ It is revealed 
that the phase boundary is smooth in accordance with 
the 3D anisotropic GL theory, and thus the system is es- 
sentially the same as the Abrikosov (or pancake) vortex 
system. A meandering phase boundary was observed in 
transport measurements by Gordeev et al. in the same 
family of materials and field regime [5J| . The difference 
might be quantitative, or may be the result of different 
experimental techniques. It is noticed that our phase 
diagram is not to scale at low fields. 

It is interesting to ask whether the phase diagram of 
interlayer Josephson vortices can be derived from a Lin- 
demann type theory. In a Lindemann theory a lattice is 
supposed to melt when the thermal average of displace- 
ment field exceeds a certain fraction of inter-vortex dis- 
tance. This picture has been useful in understanding the 
melting of pancake vortex lattice (or Bragg glass). One 
should be very careful when applying the Lindemann the- 



ory to the melting of interlayer Josephson vortex lattice, 
since thermal fluctuations are very anisotropic especially 
at strong magnetic fields above the multicritical point. 
While thermal fluctuations in the c direction are pinned 
completely below Tx , those in ab planes diverge as tem- 
perature approaches it from below. Layers of Josephson 
vortices slide easily over each other in the intermediate 
phase. A naive Lindemann theory is clearly not applica- 
ble. 

It is important to notice the differences among the 
present system and smectic liquid crystals with or with- 
out an external held. Ideal smecitc-A liquid crystals 
have liquid-like correlations in two dimensions and a 
solid like periodic modulation of the density along the 
third direction. The elastic free-energy functional for a 
smectic possesses the so-called Landau-Peierls instabil- 
ity JGJ, |62, [63|, i.e. suppression of in-plane quadratic 
first derivatives. When a magnetic or dielectric field is 
applied, the layer normal is energetically confined in a 
predetermined plane. The quadratic first derivative is 
then suppressed in only one direction. Systems in this 
group were called planar layered, and the possible phases 
and phase diagrams were addressed by Grinstein et al. 
[64| . The interlayer Josephson- vortex system in magnetic 
fields parallel to the layers is polar layered, in which the 
normal of Josephson-vortex layers is along the crystal- 
graphic c axis. There is no Landau-Peierls like instabil- 
ity left any more in the present system. The difference in 
symmetry results in the different phases, as clearly indi- 
cated by the simulated structure factors. It is interesting 
to observe the similarity between the last two layered sys- 
tems, namely both of them contain a multicritical point 
in their phase diagrams. 
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